@inbounds for i = ix:ix-1+n
    x[i] = y[i] + z[i]*a
end
